Optical spectroscopy and X-ray detections of a sample of quasars 
and AGN selected in the mid-infrared from two Spitzer wide-area 

surveys. 
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ABSTRACT 

We present optical spectroscopy of a sample of 77 luminous AGN and quasars selected on 
the basis of their mid-infrared colors. Our objects are selected from the Spitzer Extragalactic 
First Look Survey and SWIRE XMM-LSS fields, with a typical 24^m flux density of 5mJy. 
The median redshift is 0.6, with a range of ^ 0.05 — 4. Only 33% (25/77) of these objects 
are normal type-1 quasars, with no obscuration. 44% (34/77) are type-2 objects, with high- 
ionization, narrow emission lines, 14% (11/77) are dust-reddened type-1 quasars, showing broad 
lines but a dust-reddened or unusually weak quasar continuum. 9% (7/77) show no sign of 
an AGN in the optical spectrum, having either starburst spectra or spectra which could be 
of either starburst or LINER types. These latter objects are analogous to the X-ray detected 
population of AGN with weak or non-existent optical AGN emission (the "XBONGs"). 21 of our 
objects from the SWIRE field fall within moderately-deep XMM exposures. All the unobscured 
quasars, and about half the obscured quasars are detected in these exposures. This sample, 
when taken together with other samples of Spitzer selected AGN and quasars, and results from 
X-ray studies, confirms that obscured AGN dominate the AGN and quasar number counts of 
all rapidly-accreting supermassive black hole systems, at least for z < 4. This implies a high 
radiative efficiency for the black hole accretion process. 

Subject headings: quasars: general - galaxies :Seyfert - infrared:galaxies 



1. Introduction 

Quasars whose optical emission is hidden by 
dust in the optical and gas in the soft X-ray have 
always been hard to find. However, a large type- 
2 (hidden) quasar population has long been pre- 
dicted, both from the statistics of Seyfert-1 and 
2 galaxies, and from those of radio galaxies and 
radio-loud quasars. A large population of hidden 
AGN is also predicted from models of the X-ray 
background (e.g. Comastri et al. 1995; Worsley 
et al. 2005), though most models predict these 
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AGN to be of Seyfert rather than quasar lumi- 
nosity. Some success in finding dust obscured 
quasars has been obtained using selection based 
on the Two Micron AU Sky Survey (2MASS) to 
find lightly dust-rcddencd (rest frame extinctions 
Ay ~ 2) quasars, based either on 2MASS colors 
alone (Cutri et al. 2001), or on 2MASS combined 
with a radio detection and very red optical- IR col- 
ors (Glikman et al. 2005). Only recently, how- 
ever, have large samples of type-2 quasars, with 
rest-frame Ay ^ 5 been found. Zakamsaka et al. 
(2003) used the Sloan Digital Sky Survey (SDSS) 
to find objects with the high-ionization narrow- 
line spectra characteristic of type-2 AGN. Com- 
bining ISO and 2MASS data allowed Haas et al. 
(2004a) and Leipski et al. (2005) to find signif- 
icant numbers of obscured AGN, including sev- 
eral type-2 quasars. Spitzer data has been used 
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to find highly-obscured AGN, both tlirough mid- 
infrared color selection (Lacy et al. 2004 [hereafter 
Paper I], 2005a; Stern et al. 2005; Alonso-Herrero 
et al. 2005; PoUctta ct al. 2006) and joint radio- 
infrared selection (Martfnez-Sansigre et al. 2005, 
2006; Weedman et al. 2006; Donley et al. 2005). 
Typc-1 quasars selected through joint Spitzer and 
optical colors and/or morphologies have been tar- 
geted by several groups (Brown et al. 2005; Hatz- 
iminaoglou et al. 2005; Siana et al. 2006). Also, 
the mid-infrared properties of X-ray selected AGN 
have been discussed by Treister et al. (2005) and 
Barmby et al. (2005). 

Mid-infrared selection of AGN is effective be- 
cause dust heated by the AGN forms a "calorime- 
ter", whose luminosity is not as strongly affected 
by obscuration as that in most other wavelength 
bands (except the radio) (e.g. Spinoglio et al. 
1995; Meisenheimer et al. 2001: Ogle et al. 2006). 
Richards et al. (2006) show that the bolomet- 
ric correction in the mid-infrared varies little for 
different subsamples of type-1 quasars (e.g. op- 
tically or infrared luminous, optically-red or op- 
tically blue), or luminosity, and that about 15% 
of the bolometric luminosity of a quasar is emit- 
ted at mid-infrared wavelengths. Integration over 
the Richards et al. (2006) mean quasar SED as 
a whole suggests that about 40% of the emission 
from a typical quasar is reprocessed by dust in 
the infrared. This means that mid-infrared selec- 
tion of AGN is probably a good proxy for selec- 
tion on the basis of bolometric luminosity, at least 
for those objects which are not so highly obscured 
that their mid-infrared luminosity is significantly 
supressed. Much of the dust emission is thought 
to arise in the "obscuring torus" required by Uni- 
fied Schemes (e.g. Antonucci 1993), though it is 
becoming increasingly clear that host dust in the 
narrow line region of AGN also contributes sig- 
nificantly (e.g. Mason et al. 2006). AGN are an 
important global contributor to the mid-infrared 
luminosity density. Brand et al. (2006) show that 
the luminosity density of mid-IR emission in sam- 
ples of MIPS 24/xm sources selected at mJy levels 
is dominated by AGN, only at submJy levels do 
star forming galaxies begin to dominate. 

In this paper, we focus on optical spectroscopic 
follow-up of AGN candidates selected according 
to the mid-infrared color criteria presented in Pa- 
per I. By selecting relatively bright objects from 



wide-area Spitzer surveys we arc able to find many 
obscured AGN with mid-infrared luminosities typ- 
ical of quasars at moderate redshifts. These ob- 
jects provide a good complement to the mostly 
lower-luminosity or higher redshift AGN found in 
surveys such as the Great Observatories Origins 
Deep Survey (GOODS) and the Chandra Deep 
fields. Preliminary results from the spectroscopic 
follow-up of the sample presented in Paper 1 were 
given in Lacy et al. (2005a). In this paper, we de- 
scribe the selection of samples of mid-infrared se- 
lected AGN candidates flux-limited at 24/Ltm from 
both the Extragalactic First Look Survey (XFLS) 
and the Spitzer Wide-area Infrared Extragalactic 
Survey (SWIRE) XMM-LSS field (Lonsdale et al. 
2005), and their spectrosopic follow-up observa- 
tions. We classify the optical spectra according to 
the type of AGN, and use them to place a lower 
bound on the fraction of obscured quasars and lu- 
minous Seyfert galaxies. We also investigate the 
X-ray fluxes of those objects which fall on archival 
observations made with the X-ray Multiple Mirror 
telescope [XMM- Newton). A more detailed anal- 
ysis of the Spectral Energy Distributions (SEDs) 
and X-ray properties of the AGN will be carried 
out in a later paper (Petric et al. 2006). 

We assume a cosmology with f^M — 0.3, Oa = 
0.7 and Hq = TOkms-^Mpc-^ 

2. The Spitzer XFLS and SWIRE AGN 
samples 

We selected samples of AGN by mid-infrared 

color from the XFLS and SWIRE XMM fields. 
The basis of our color selection technique is de- 
scribed in Sajina, Lacy &: Scott (2005) and Pa- 
per I, though the samples in this paper are fiux- 
limited at 24/im rather than at 8/xm as was the 
case for the sample of Paper 1. To obtain the 
XFLS sample, we began by matching the XFLS 
24/xm catalog of Fadda et al. (2006) to the In- 
frared Array Camera (IRAC) 4-band catalog of 
Lacy et al. (2005). All MIPS 24/^m detections 
with detections in all four IRAC bands were then 
plotted on the color-color diagram of Figure 1. 
For the SWIRE XMM sample, we used the band- 
merged SWIRE IRAC/MIPS 24/im catalog avail- 
able from the Spitzer popular products website. 
The color cut used to select AGN candidates was 
lg(55.8/^3.6) > -0.1 and lg(58.o/54.5 > -0.2 
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and lg(58.o/54.5) < 0.81g(55.8/^3.6) + 0.5. The 
lg(5'5.8/'S'3.6) cut was moved by +0.1 compared to 
Paper I. The exact position of this cut is not criti- 
cal, and it was fch that shifting this would help re- 
move non-AGN contaminants, particularly as the 
final 5.8/xm fluxes in the XFLS catalogue were in- 
creased by « 10% in the final version of the XFLS 
catalog relative to the catalog used in Paper I. We 
discuss the exact placement of the color cuts with 
respect to the completeness of the AGN selection 
in section 5.3. Candidates were selected down to 
limiting 24/Ltm fluxes of 4.4mJy in the XFLS and 
6.6mJy in the SWIRE XMM field. 

The 24/im flux limited sample should be an im- 
provement on 8/xm selection by making us more 
sensitive to highly-reddened objects, and high red- 
shift objects, where the strong fc-correction on 
the mid-infrared dust emission means that they 
drop out of flux-limited samples selected at shorter 
wavelengths. The disadvantage of this selection is 
that we will be missing objects at z ~ 1.5 which 
have strong silicate absorption in the 24/im band, 
however, as our median redshift is « 0.6 only a 
few objects will be missed due to this eff'ect. The 
XFLS sample is complete for objects falling in the 
IRAC "AGN wedge" color selection to a flux limit 
of 4.6mJy at 24/zm, and the SWIRE sample to a 
flux limit of 6.6mJy (note that all MIPS sources 
down to these flux limits have 4-band IRAC detec- 
tions). Approximately 20% of MIPS 24/im sources 
at these flux levels fall into our AGN-dominated 
part of the color-color plot. 

A few stars and nearby galaxies have been re- 
moved from the samples. Most stars were only 
present in the initial selection due to photomet- 
ric errors related to saturation in the shorter IRAC 
bands, though one star (SSTXFLS J172432. 8+592646, 
sec Paper I) has the correct colors to fall in the 
sample, presumably due to a dust disk. Very 
low redshift {z < 0.05) galaxies are also possi- 
ble contaminants due to the presence of 6.2/im 
PAH emission in the 5.8/im band of IRAC (Sajina, 
Lacy & Scott 2005). The interacting galaxy Arp 
54 appeared in the selection of the SWIRE AGN 
sample, but has been removed for this reason. 

Tables 1 and 2 contain the observing logs for 
our samples. Most data were taken with the Cos- 
mic instrument on the Palomar 200-inch. Some 
spectra were obtained with the ISIS instrument 
at the 4.2m WiUiam Herschel Telescope (WHT) 



on La Palma, with ESI at the Keck II tele- 
scope, or with the Kast spectrograph at the 3m 
Shane Telescope at Lick Observatory. SSTXFLS 
J172248. 9+583256, for which we failed to obtain 
a redshift in the optical was observed in the near- 
infrared using the low dispersion mode of the 
Spcx instrument at the Infrared Telescope Facility 
(IRTF) (Rayner et al. 2003). Data were analysed 
in IRAF using standard routines for flat flelding, 
wavelength calibration and sky subtraction. No 
attempt has been made to correct for redden- 
ing. An extinction correction has been applied 
based on airmass, but no attempt has been made 
to correct the telluric absorption features. The 
spectra are approximately photometrically cali- 
brated, but variable seeing and cirrus cloud means 
that we expect our spectrophotometric accuracy 
to be no better than 50%. In the XFLS field, 
most quasars have spectra from the Sloan Digital 
Sky Survey (SDSS), which have been used where 
they exist. A few of the quasars in the SWIRE 
field have redshifts from the literature. Several 
of our XFLS sources have spectra in the Multiple 
Mirror Telescope (MMT) spectroscopic survey of 
Papovich et al (2006; hereafter P06); including 
one object missed from our original selection but 
which is within the flux-limited sample (SSTXFLS 
J171454.4+584948). The spectra of aU the ob- 
jects with the exception of the SWIRE quasars 
with redshifts from the literature are shown in 
Figures 2 and 3. Objects with spectral classifi- 
cations from both this paper and Papovich et al. 
(2006) are discussed in Appendix A. For complete- 
ness. Appendix B contains details of the results of 
spectroscopy on the objects in the 8/um sample of 
Paper I which are too faint at 24/im to fall in the 
XFLS sample. 

3. Results 

3.1. Optical spectral classification of the 
AGN candidates 

We have used a combination of standard classi- 
fication schemes to classify our AGN candidates, 
principally because the large redshift and luminos- 
ity range of our sample makes any single scheme 
unworkable. We divide our objects into four main 
classes. Normal quasars/Scyfcrt-1 galaxies are as- 
signed type 1. The obscured quasars are subdi- 
vided into three classes. Reddened quasars, show- 
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Table 1: Observing log for the XFLS AGN sample 
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95 
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90 
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1.5 


1x90076 
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2005 Juno 02 
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1.5 


1.15 
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SSTXFLS 
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Palomar 200/Cosmic 


2x600 
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90 
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90 
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1.13 
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1.18 


45 
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Shane /Kast 


Ix 1200 


1.5 


1.24 


73 


SSTXFLS 


J171530. 7+600216 


2005 June 01 


Palomar 200/Cosmic 


1x600 


1.5 


1.15 


140 


S 

SSTXFLS 


,1171650.6+595752 


2005 Juno 02 


Palomar 200/Cosmic 


1x600 


1.5 


1.23 
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90 
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90 
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1x600 
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1.57 
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Palomar 200/Cosmic 
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1.48 


90 
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J172228.1+601526 
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Palomar 200/Cosmic 


1x1200 


1.5 


1.24 


45 


SSTXFLS 
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2005 Nov 01 


Palomar 200/Cosmic 


2x1200 


1.5 


1.55 
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SSTXFLS 


J172248. 9+583256 


2005 June 02 


Palomar 200/Cosmic 


2x1200 


1.5 


1.12 


45 


SSTXFLS 
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0.5 


1.48 
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2004 August 19 
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90 
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1.34 
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Table 2: Observing log for the SWIRE XMM-LSS sample 
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7 
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42 


SWIRE2 


J021939, 


.08-051133. 
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29 
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90 
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01 


Palomar 


200/Cosmic 


1x600 


1, 
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47 
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47 
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Fig. 2. — Optical spectra of the XFLS AGN candidates. Dotted lines indicate the positions of redshifted 

emission lines (not all of which are necessarily detected in a given spectrum). The emission lines indicated are 
(from blue to red): Lya, Nvl240, Civ 1549, Heiil640, Ciii]1909, [Neiv]2424, Mgii2798, [Nev]3426, [Oii]3727, 
[Neiii]3869, H7, H/3, [Oiii]4959, [Oiii]5007, [Nii]6548, Ha and [Nii]6584. The spectra are presented in order 
of Right Ascension. 
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2. — Optical spectra of the XFLS AGN candidates, continued 
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2. — Optical spectra of the XFLS AGN candidates, continued 
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2. — Optical spectra of the XFLS AGN candidates, continued 
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Fig. 3. — Optical spectra of the SWIRE AGN candidates. The position of emission lines are indicated 
described in the caption to Figure 2. The spectra are presented in order of Right Ascension. 
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3. — Optical spectra of the SWIRE AGN candidates, continued 
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3. — Optical spectra of the SWIRE AGN candidates, continued 
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Fig. 4. — Two diagnostic (BPT) diagrams for our lower redshift emission line galaxies that have Ha in the 
optical window. The dividing lines are from Kewley et al. (2006). Most of our objects plot amongst the 
AGN, there are no LINERs in the sample. 



' Type- 


1 


A Type- 


IR 


• Type- 


2 


Sbt 4: 


Sbt/L 



lg(S5.„/S3.e) 

Fig. 1.— IRAC color-color plot for the XFLS field 
for objects detected in MIPS 24/im and all four 
IRAC bands. The region within the dashed line is 
that used to select AGN. Objects in the AGN sam- 
ples (both XFLS and SWIRE) are plotted as col- 
ored symbols. Object classifications are discussed 
in Section 3. 



ing broad lines in their optical spectra, but with 
a faint, presumably dust-reddened quasar contin- 
uum and/or clear host galaxy emission with no 
sign of a blue quasar continuum above the host 
galaxy light are classed type-lR. Objects with 
high-ionization narrow emission lines which can- 
not be produced in a pure starburst scenario are 
classed as type-2. Objects with starburst emis- 
sion line ratios in emission line diagnostic dia- 
grams (Baldwin, Phillips & Terlevich 1983; here- 
after BPT), and/or low ionization spectra with 
strong Balmer lines are classed as starbursts (Sbt). 
Where we have a low-ionization spectrum with 
too few lines to plot on a diagnostic diagram, 
or are low ionization and have Ha beyond the 
range of the optical spectrum they are classed as 
starburst/LINER (Sbt/L). Details of the bases of 
our optical classifications are described below. It 
should be emphasized that these classes are based 
on the observed-bame optical spectra, the rest- 
frame optical spectra may give different classifi- 
cations. The classifications are given in the final 
columns of Tables 3 and 4, together with their 
basis. Of our sample of 77 objects, we find 25 ob- 
jects with type-1 optical spectra, 11 with type-lR, 
34 with type-2, two with unambiguous starburst 
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spectra and five with spectra which could be from 
either starbursts or LINERs. 

3.2. Type-1 objects 

Objects with Unewidths of at least one permit- 
ted emission line > 3000km~^ are classified as 
type- Is. These objects have "BL" in the classifica- 
tion basis columns in Tables 3 and 4. This is suffi- 
cient to distinguish most of the luminous quasars. 
However, as our sample spans a wide luminos- 
ity range, several of the lower luminosity, lower 
redshift type-1 objects have narrower permitted 
lines (including a couple of candidate narrow-line 
Seyfert-1 galaxies). We therefore use the presence 
of strong permitted lines (H/3/[Oiii]5007> 0.5) 
and Fell emission in the 4000-5500A range to clas- 
sify them as type-1. These objects have "Fell" in 
the basis columns. The division into normal type- 
1 or reddened, type-lR, quasars is also made. This 
is decided on one of two criteria. If the optical 
spectrum is quasar-like, with featureless contin- 
uum and broad emission lines, but the continuum 
appears significantly reddened {E(B — V) ^ 0.3) 
compared to that of a normal quasar (see section 
3.5), then the object is classified at type-lR. Also, 
some low-z objects were classed as type-lR if the 
optical spectrum contains broad lines, but host 
galaxy features are visible and the object is dom- 
inated by host galaxy light in the blue. 

3.3. Type-2 objects 

Objects without broad lines or Fell emission, 
but with high-ionization narrow lines are classi- 
fied as follows. Objects with [Nev] detections with 
rest-frame equivalent width > 5A are classed as 
type-2 with the basis column containing "[Nev]" 
and the measured rest-frame equivalent width in 
Angstroms (for comparison the [Nev] 3426 emis- 
sion line in the SDSS composite quasar spec- 
trum of Vanden Berk et al. (2001) has an EW 
of « 0.6A). For low redshift objects where Ha 
was included in the spectra, we have used BPT 
diagrams to classify the AGN. Figure 4 shows 
the emission line ratios of our objects plotted on 
two diagnostic diagrams using the divisions be- 
tween star-forming, Seyfert and LINER galax- 
ies proposed by Kewley et al. (2006). Most ob- 
jects plot with the Seyferts, though two objects 
plot close to the AGN/starburst division line in 



both plots, and are presumably composite star- 
burst/AGN objects (SSTXFLS 171650.6-1-595752; 
SWIRE2 J022225.86-050015.1). We class these as 
type-2 as they do have evidence of AGN activity, 
even if this is not dominating their optical emission 
line ratios. For some objects with neither [Nev] 
nor Ha/ [Nil] detected in their spectra we have re- 
lied on the ionization parameter as an approximate 
AGN indicator and used the [Oiii] 5007/ [On] 3727 
ratio, classing objects with [Olll] 5007/ [Oil] 3727 
> 1 as questionable type-2s. 

There are also a few high redshift objects 
having high-ionization UV permitted lines with 
widths ~ 2000kms~^. Linewidths of this order are 
seen in the UV spectrum of NGC1068 (Grimes et 
al. 1999). Given this, and their high luminosities, 
we therefore anticipate that the true broad line 
widths will be much higher, so we are probably 
seeing only the narrow-line region in our spectra. 
However, given their high redshifts {z > 1.5), only 
a very small amount of extinction is required to 
hide the broad-line region in the rest-frame UV. 
For the three highest redshift objects, at 2.78, 2.99 
and 4.27, the rest-frame optical to near-infrared 
light is sampled in the IRAC bands, and this must 
almost certainly be dominated by reddened quasar 
light for these objects to exceed our IRAC flux lim- 
its. It is therefore unclear whether these objects 
should be classed as type-2s (showing no strong 
broad emission lines or quasar continuum in the 
rest-frame optical), or as class IRs. We thus as- 
sign z > 1.5 objects with narrow UV lines class 
2HZ (and group them with the type-2 objects for 
plotting and statistical purposes), pending better 
information on their rest-frame optical emission 
line properties and SEDs from near infrared imag- 
ing and spectroscopy. In all of these objec;ts, Lya 
is typically faint compared to the high ionization 
UV lines, probably due to absorption by neutral 
hydrogen (e.g. van Ojik et al. 1994). 

3.4. Starbursts and LINERs 

We are left with seven objects with low- 
ionization lines, showing no signs of obvious AGN 
activity. Such objects have been classified as 
"type-3" AGN by Leipski et al. (2005). They 
are all relatively low redshift (z < 0.6). All show 
[Oll]3727 in emission, and Balmer absorption fea- 
tures. One also shows strong narrow H/3 and H7 
emission. The two objects we have been able to 
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positively identify as starburst in the optical have 
been classified as such. The remainder are clas- 
sified as starbursts/LINERs. Note that we fail 
to find any unambiguous LINERs in our sam- 
ple. (SSTXFLS 171454.4+584948 plots below the 
Seyferts in the N[ii]/Ha vs [Oiii]/H/3 diagram, 
in the traditional LINER region [e.g. KaufFmann 
et al. 2003] but well within the starbursts in the 
[Oi]/Ha vs [Oiii]/H/3 diagram which gives a bet- 
ter discrimination of LINERs according to Kew- 
ley et al. [2006], so we have classified it as a star- 
burst.) This may be due to LINERs lacking a sub- 
stantial hot torus bright enough to dominate the 
mid-infrared over stellar emission (e.g. Ogle et al. 
2006), but better spectra of our starburst/LINER 
objects are needed to confirm this. 

For the most part our classifications are consis- 
tent with those of Papovich et al. (2006) , though 
there are some objcicts which are differently clas- 
sified. These are discussed in Appendix A. 

3.5. Optical to mid-infrared color 

In Figure 5 we show the observed B - [8.0] color 
plotted against redshift for the different spectral 
classifications. For the XFLS sample, the R mag- 
nitudes are from the survey of Fadda et al. (2004) , 
those for the SWIRE sample are from the Cam- 
bridge Automatic Plate Measuring (APM) scans 
of the UK Schmidt Telescope plates (though in 
a few cases of faint or blended sources the mag- 
nitudes have been estimated by eye from digi- 
tized sky survey plates). Apart from a couple of 
low redshift AGN, which probably have their op- 
tical emission dominated by their host galaxies, 
all the normal type-1 objects trace the expected 
color locus for type-1 quasars very well, particu- 
larly the "IR luminous" BED of Richards et al. 
(2006), suggesting that the mid- infrared selection 
is not picking a very unusual population of type- 
1 quasars, but rather is finding quasars with the 
properties that wc would expect from a sample se- 
lected in the mid-infrared. After the host galaxy 
contribution to the i?-band becomes negligible at 
z ^ 0.4, the normal and obscured (type IR, 2 and 
3) quasars separate, with very few objects hav- 
ing moderately-red colors. This relatively clean 
separation shows that our optical spectral classifi- 
cations are indeed likely to be correct. 

The z > 1.5 type-2 objects (type 2HZ in Ta- 
ble 3), which have intermediate UV linewidths 



and whose precise optical classification is therefore 
unclear, are obviously redder than the expected 
mid-infrared to optical color of normal quasars at 
these redshifts, suggesting that they are not sim- 
ply quasars with unusually narrow lines, but do 
indeed have most of their nuclear emission extin- 
guished by dust. 

Our division between type-1 and type-lR 
quasars also appears to be justified on the basis 
of this color (at least for objects at z s; 0.4 where 
the host galaxy does not affect the i?-magnitude), 
with the type-lR quasars having rest-frame red- 
denings E{B - F) ~ 0.3 - 0.7. 

4. X-ray detections 

There is some overlap between objects in the 
SWIRE-XMM fields and archival XMM-Newton 
data. 21 of 34 sources in the SWIRE-XMM sample 
fall on moderately deep (5-50ks) XMM exposures. 
These are sufficiently deep to detect, or to set 
useful limits on, the X-ray flux from our sources. 
These data appear to have been taken for sev- 
eral different surveys, including the XMM large- 
scale structure survey (e.g. Chiapetti et al. 2005) 
and the the Subaru XMM-Newton Deep Survey 
(SXDS; Watson et al. 2004). We have obtained 
approximate fluxes for these objects in order to in- 
vestigate the optical-to-X-ray luminosity ratios of 
our sources (Table 5). This allows us to place our 
mid-infrared selected obscured AGN into the con- 
text of X-ray samples of similar objects. Counts 
in the 0.2-12keV range from the PN detector of 
the European Photon Imaging Camera (EPIC), 
as listed in the IXMM serendipitous catalog (e.g. 
Michel et al. 2004), or the pre-release of the 2XMM 
catalogue (The Second XMM-Newton Serendipi- 
tous Source Pre-release Catalogue, XMM-Newton 
Survey Science Centre) have been used where 
available, otherwise counts have been taken from 
the source lists available from the XMM archive 
(the pipeline products [PPS] lists). Limits are 
estimated from the fluxes of the faintest believ- 
able detections on each image, and correspond to 
Ri 3(7. Fluxes from the IXMM or 2XMM catalogs 
are from the PN countrate in IXMM, and from the 
average over all detectors in 2XMM. Where only 
counts or a limit are available fluxes are derived 
from the count rate using the PIMMS website 
(http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html). 
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Table 3: Bright AGN in the XFLS region 



Name 


S24.o/mJy^ 


Redshift 


Class 


Basis 


SDSS J171033. 1+584456 


6.1 


0.281 


IR 


Fcii/BL 


SSTXFLS J171115. 2+594906 


9.4 


0.587 


Sbt 


strong Balmcr lines 


SDSS J171117.6+584123 


5.8 


0.617 


1 


BL 


SDSS J171124.2+593121 


5.5 


1.491 


1 


BL (LoBAL?) 


SSTXFLS J171147.4+585839 


4.8 


0.800 


2 


[Nev]13 


SDSS J171207.4+584754 


13.4 


0.269 


1 


Fell 


SSTXFLS J171233. 4+583610 


5.1 


1.663 


IR 


BL 


SSTXFLS J171302. 3+593610 


11.8 


0.668 


2 


[Ncv]33 


SSTXFLS J171313. 9+603146 


10.5 


0.124 


2 


BPT 


SSTXLFS J171324.2+585549 


4.9 


0.609 


2 


[Nev]12 


SSTXFLS J171325. 1+590531 


9.5 


0.126 


2 


BPT 


SSTXFLS J171331. 5+585804 


5.9 


0.435 


Sbt/L 


[Oin]/[On]<l 


SDSS J171334.0+595028 


5.4 


0.615 


1 


BL 


SSTXFLS J171335.1+584756 


23.7 


0.133 


IR 


Fell 


SSTXFLS J171345.5+600730 


4.9 


2.78 


2HZ 


narrow UV lines 


SDSS J171352. 4+584201 


23.9 


0.521 


1 


BL/Fcii 


SSTXFLS J171419. 9+602724 


5.6 


2.99 


2HZ 


narrow UV lines 


SSTXFLS J171430. 7+584225 


8.2 


0.561 


IR 


BLt 


SSTSFLS J171454.4+584948* 


4.9 


0.253 


Sbt 


BPT 


SSTXLFS J171513.8+594638 


5.0 


0.248 


IR 


BL (radio-loud) 


SSTXFLS J171530. 7+600216 


11.6 


0.420 


2 


[Nev]26 


SSTXFLS J171650. 6+595752 


6.7 


0.182 


2 


BPT (composite) 


SSTXFLS ,7171708.6+591341 


5.4 


0.646 


2 


[Nev]37 


SDSS J171732. 9+594747 


4.6 


1.058 


1 


BL 


SDSS J171736. 9+593011 


6.4 


0.599 


1 


BL 


SDSS J171747.4+593258 


5.4 


0.248 


1 


BL 


SSTXFLS J171750.7+584745 


5.1 


2.55? 


IR? 


BL? 


SSTXFLS J171754.6+600913 


8.9 


4.27 


2HZ 


narrow UV lines 


SSTXFLS J171831. 5+595317 


8.3 


0.700 


2 


[Nev]26 


SSTXFLS J171839.6+593359 


10.9 


0.383 


1 


BL (radio loud) 


SDSS J171902. 2+593715 


26.9 


0.178 


1 


BL/Feii 


SSTXFLS J171913. 5+584508 


8.9 


0.318 


2 


[Nev]12, BPT 


SSTXFLS J172044.8+582923 


5.1 


1.697 


2HZ 


narrow UV lines 


SSTXFLS J172050.4+591511 


9.4 


1.44 


IR 


BL 


SSTXFLS J172123. 1+601214 


13.3 


0.325 


2 


[Nev]33/BPT 


SSTXFLS J172219. 5+594506 


7.8 


0.271 


2 


BPT 


SSTXFLS J172228. 1+601526 


7.1 


0.741 


2 


[Nev]15 


SDSS J172238. 7+585107 


6.7 


1.617 


1 


BL 


SSTXFLS ,7172245.0+590328 


4.7 


0.797 


2 


[Nev]43 


SSTXFLS J172248. 9+583256 


5.6 


0.84 


IR 


BL 


SSTXFLS J172328.4+592947 


8.1 


1.34? 


IR 


BL 


SSTXFLS J172542.3+595317 


9.5 


0.437 


Sbt/L 


low ionization 


SDSS J172619.7+601600 


6.6 


0.925 


1 


BL 



* Spectrum from Papovich et al. (2006). 

^ classification based on Papovich et al. spectrum, see Appendix A. ■'■ Optical/near-IR redshift is uncertain, but a Spitzer mid-IR 
spectrum shows silicate absorption at z = 0.84 (Lacy et al. 2006, in prepraration). 
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Tabic 4: Bright AGN in the SWIRE-XMM region 



Name 




S'24.o/mJy 


Rcdsliifl 


Class 


Basis 


SWIRE2 


J021638. 21-042250. 8 


14.6 


0.304 


2 


BPT 


SWIRE2 


J021640.72-044405.1 


14.8 


0.870 


1 


Becker et al. (2001) 


SWIRE2 


J021657.77-032459.7 


24.2 


0.137 


1 


BL 


SWIRE2 


J021729. 06-041937. 8 


8.9 


1.146 


IR 


BL 


SWIRE2 


J021749.00-052306.9 


8.1 


0.987? 


Sbt/L? 


low ionization? 


SWIRE2 


,7021759.91-052056.1 


8.8 


0.535 


2 


[Ncv]26 


SWIRE2 


J021808. 22-045845. 3 


9.2 


0.712 


1 


Caccianiga et al. (2004) 


SWIRE2 


J021809. 45-045945. 9 


7.8 


1.094 


2 


[Nev]90 


SWIRE2 


J021830.57-045622.9 


8.5 


1.401 


1 


Sharp et al. (2002) 


SWIRE2 


J021859.74-040237.2 


16.0 


0.203 


2 


BPT 


SWIRE2 


,1021909.60-052512.9 


25.6 


0.099 


2 


BPT 


SWIRE2 


,7021938.70-032508.2 


6.6 


0.435 


1 


Fell 


SWIRE2 


,7021939.08-051133.8 


32.9 


0.150 


2 


[Nev]9 


SWIRE2 


,7022005.93-031545.7 


6.6 


1.560 


1 


BL 


SWIRE2 


J022012. 21-034111. 8 


6.8 


0.166 


1 


BL 


SWIRE2 


J022039.48-030820.3 


10.6 


0.451 


1 


BL 


SWIRE2 


J022133.82-054842.8 


6.9 


0.502 


1 


BL 


SWIRE2 


J022211. 57-051308.1 


6.7 


0.841 


2 


[Nev]30 


SWIRE2 


J022225. 86-050015.1 


8.3 


0.149 


2 


BPT (composite) 


SWIRE2 


J022255.87-051351.7 


8.3 


0.846 


1 


BL 


SWIRE2 


,7022301.97-052335.8 


6.8 


0.708? 


Sbt/L? 


low ionization? 


SWIRE2 


J022306. 74-050529.1 


16.0 


0.330 


2 


BPT 


SWIRE2 


J022310.42-055102.7 


13.3 


0.506 


2 


[Nev]15 


SWIRE2 


J022348.96-025426.0 


7.6 


0.355 


3 


low ionization 


SWIRE2 


J022356.49-025431.1 


10.6 


0.451 


2 


[Nev]16 


SWIRE2 


,7022422.27-031054.7 


7.8 


1.226 


1 


BL 


SWIRE2 


,7022431.58-052818.8 


9.4 


2.068 


1 


BL 


SWIRE2 


,7022438.97-042706.3 


6.6 


0.252 


2 


BPT 


SWIRE2 


,7022455.77-031153.6 


7.6 


0.790 


2? 


[Oiii]/[On]>l 


SWIRE2 


J022508. 33-053917.7 


9.6 


0.293 


2 


BPT 


SWIRE2 


J022542. 02-042441.1 


10.9 


0.580 


2 


[Nev]23 


SWIRE2 


J022600.01-035954.5 


9.0 


1.095? 


IR 


BL 


SWIRE2 


J022612. 67-040319. 4 


7.7 


0.330 


2 


BPT 


SWIRE2 


J022700.77-042020.6 


27.1 


0.053 


1 


Fen 
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assuming an input spectrum with photon index of 
two and a column of lO^^cm"^. 

A few points are clear from even this prelim- 
inary study. All the type-1 quasars are detected 
in the X-ray (albeit with about an order of mag- 
nitude of scatter in their X-ray to infrared ra- 
tios), but only about half of the optically-obscured 
quasars (Figure 6). The obscured quasars do, how- 
ever, show a wide range in X-ray to infrared flux 
ratios, ranging from little or no absorption to X- 
ray fluxes absorbed by a factor of ten or more 
(averaged over the 0.2-12keV band) relative to 
the type-Is. Interestingly, one of the two star- 
burst/LINER objects in these X-ray fields is de- 
tected, though the other is not. 

5. Discussion 

5.1. Selection effects and redshift biasses 

Every technique for selecting AGN is affected 
by selection biases, and this one is no exception. 
The color selection means that objects whose ob- 
served mid-infrared colors are not dominated by 
thermal emission from the AGN will be missing 
from the sample. This can happen in a variety of 
ways. First, the AGN may simply be of such low 
luminosity in the infrared that the starlight from 
the host galaxy dominates the mid-infrared emis- 
sion, either because the AGN is intrinsically weak, 
or because the covering factor of the hot dust is 
low, perhaps because of an unusual torus geome- 
try. This bias gets worse with increasing redshift 
as the stellar bump in the rest-frame near-infrared 
is redshifted into the mid-infrared bands, and the 
fc-correction on the thermal emission of the dust 
becomes large. Second, some AGN will simply 
be so highly-obscured that the AGN emission is 
not seen at observed wavelengths < 8/im. Third, 
some AGN may also be so highly obscured that 
they drop out of the IRAC catalogs at short wave- 
lengths. 

The positions of the different classes of objects 
within the "AGN wedge" on the IRAC color-color 
plot of Figure 1, shown in detail in Figure 7, can 
give us some idea of the particular selection effects 
for each class of object. The type-1 quasars all lie 
in a fairly tight group in the center of the wedge, 
so few are likely to be missing (see also Paper 1). 
Reddened type-Is show more dispersion, includ- 
ing two with very red IRAC colors. The type- 
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Fig. 5. — Optical to infrared color for our quasar 
samples. Type-3 objects which lack redshifts have 
been assigned a redshift of 1.50. The color loci of 
two of the quasar SEDs of Richards et al. (2006) 
are plotted as lines. The blue solid line is the mean 
quasar SED and the green long-dashed line the 
"IR luminous" SED. The red dashed lines indicate 
the colors of the mean quasar SED reddened by 
the Small Magellenic Cloud extinction law of Pei 
(1992). The lower line corresponds to a rest-frame 
reddening of E(B — 1/) = 0.32 and the upper to 
0.65. The dispersion in optical to mid-IR color of 
the quasars of Richards et al. (2006) is « ±0.3 
mag. 
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Fig. 6. — Observed X-ray (0.2-12 keV) flux versus 
observed S.O/Ltm infrared flux ratios for objects in 
archival XMM fields. The blue solid line is the 
mean quasar SED of Richards et al. extrapolated 
from the soft X-ray assuming a photon index of 
1.8, the green dashed line the same for the "IR 
luminous" SED. 
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Fig. 7. — Enlarged view of the "AGN wedge" of 
Figure 1. Large symbols represent z > 0.8 objects, 
small symbols those at z < 0.8. Open squares 
represent the objects from the sample of Table 6 
which have AGN-like optical emission line ratios 
but non-AGN IRAC colors. 
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Fig. 8.— [Oi]/Ha vs [Oiii]/H/3 BPT diagram for 
the objects of Table 6 (the 50 brightest 24/im 
sources, regardless of IRAC colors) with [Oili] de- 
tections. Objects falling outside of the IRAC AGN 
color selection are shown as open symbols, those 
within it as closed. 



19 



2s are also more dispersed, with several (mostly 
low redshift objects) trending towards the diago- 
nal divide in the color-color plot, presumably due 
to PAH emission, and several with very red colors. 
The starbursts and starburst/LINERs are surpris- 
ingly close to the center of the wedge, though one 
does have extremely red IRAC colors. From Fig- 
ure 7 we can deduce that some (probably mostly 
lower-2) type-2 objects are probably missing due 
to plotting with the PAH-dominated galaxies, and 
that all classes of obscured quasars are probably 
missing their reddest representatives. There is a 
weak trend whereby the higher redshift objects 
tend to be located towards redder colors in the 
wedge, though at the highest redshifts, where the 
IRAC colors are sampling the rest-frame optical 
and near-infrared, the colors become bluer again. 

As these selection effects have yet to be modeled 
in depth, we will defer further discussion to a later 
paper, where we will attempt to model the emis- 
sion from our objects in detail and apportion lu- 
minosity between AGN and starburst components 
(Petric et al. 2006). However, we point out here 
that our fraction of obscured quasars is necessarily 
a lower limit, and this limit becomes increasingly 
poor with increasing redshift. 

5.2. Objects with low ionization optical 

spectra 

Are there any high luminosity AGN in which 
both the broad and narrow-line regions are ob- 
scured, or otherwise missing? This study has 
found several objects which qualify as "type-3" 
quasars (following Leipski et al. [2005]), namely 
galaxies with emission line properties more typical 
of LINER or starburst galaxies than narrow-line 
AGN. Leipski et al. find 57% of their IS0/2MASS 
selected AGN have these optical spectra. Af artinez- 
Sansigre et al. (2005, 2006) also deduce the 
presence of a significant population of obscured 
narrow-line AGN from their mid-infrared/radio 
selected objects which show no emission lines in 
the rest-frame UV. 

Our study has a lower fraction. Including 
both objects with starburst or ambiguous star- 
burst/LINER spectra we find a "type-3" fraction 
of w 9%, which is of course an upper limit as, apart 
from one case which is detected in the X-ray at the 
level expected for an AGN (SWIRE2 J021749.00- 
052306.9), we cannot prove the presence of AGN 



in these objects. Thus, despite the inferred pres- 
ence of hot dust, infrared spectroscopy or X-ray 
detection (at a level greater than expected for a 
starburst) is required to prove the AGN nature 
of these candidate type-3 objects. If they are red- 
dened AGN, the reddening may be in the host (e.g. 
Martmez-Sansigre et al. 2005; Rigby et al. 2006), 
though we cannot eliminate the possibility that 
nuclear dust is obscuring the narrow-line region, 
or is blocking the line of sight to the narrow line 
region from the nucleus. 

X-ray surveys too find a population of AGN 
which lack opical evidence for AGN activity, the 
X-ray bright, optically- normal galaxies (XBONGs) 
(e.g. Comastri et al. 2002). A number of processes 
have been suggested for the lack of optical AGN 
activity in these objects, including host galaxy 
obscuration (Rigby et al. 2006), radiatively inef- 
ficient accretion (Yuan & Narayan 2004) or dilu- 
tion of weak emission lines by host galaxy light 
(Moran, Filippenko & Chornock 2002; Eckhart et 
al. 2006). The discovery of [Nv]3426 emission in a 
stacked spectrum of XBONGS made by Eckart et 
al. (2006) favors dilution and/or extinction over 
a fundamental difference in accretion mode. The 
same may well be true if our type-3 quasars re- 
ally are AGN, particularly any that are composite 
starburst/ AGN systems with strong star forma- 
tion activity in the host galaxies. In particular, 
objects which are classified on the basis of line 
ratios in the BPT diagram may have their AGN 
line emission overwhelmed by the starburst com- 
ponent, which is likely to dominate the overall 
energetics of the galaxy. 

Dilution of AGN signatures by host galaxy 
emission may help to explain why our type-3 frac- 
tion is lower than those of Leipski et al. , whose ob- 
jects are mostly of lower rcxlshift and luminosity, 
where host galaxy contamination is more of an is- 
sue. Our type-3 fraction may be smaller than that 
of Martiinez-Sansigre et al. (2005,2006) due to our 
order of magnitude brighter flux limit at 24/xm al- 
lowing us to find more objects with intrinsically 
faint or lightly-reddened {E{B — V) ^ 0.1 — 1) nar- 
row lines. Certainly it is apparent from our spec- 
troscopic campaign that even objects with 24yum 
fluxes of several mJy can require 1.5hrs of integra- 
tion on a 5m telescope to see emission lines. 
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Tabic 5: XMM-Newton detections of AGN in the SWIRE field 

Object name z Opt Obs ID Exposure Camera counts source 'S'o,2-i2kcv 

class name time ( 1 ~ 

ergs~"'"cm~^) 



SWIR.E2 J021749. 00-052306. 9 


0.987? 


Sbt/L?? 


0112370701 


47648 


PN 


20050 


2XMM 


178 


SWIR.E2 J021759. 91-052056.1 


0.535 


2 


0112370801 


47330 


PN 


< 50 


PPS 


< 0.4 


SWIRE2 ,7021808.22-045845.3 


0.712 


1 


0112370101 


57631 


PN 


14500 


IXMM 


104 


SWIRE2 J021809. 45-045945. 9 


1.094 


2 


0112370101 


57631 


PN 


100 


PPS 


0.75 


SWIRE2 J021830.57-045622.9 


1.401 


1 


0112370101 


57631 


PN 


3100 


IXMM 


22 


SWIRE2 ,7021909.60-052512.9 


0.099 


2 


0112370801 


47330 


PN 


300 


2XMM 


2.6 


SW1RE2 ,7021939.08-051133.8 


0.150 


2 


0112370801 


47330 


PN 


601 


IXMM 


5.3 


SW1RE2 J022012. 21-034111. 8 


0.166 


1 


0037982601 


11592 


PN 


1740 


2XMM 


63 


SW7RE2 J022039. 48-030820. 3 


0.451 


1 


0037982701 


14892 


PN 


1100 


2XMM 


31 


SW1RE2 ,7022211.57-051308.1 


0.841 


2 


0111110501 


20056 


PN 


;S 20 


PPS 


S 0.4 


SWIRE2 J022225. 86-050015.1 


0.149 


2 


0111110501 


20056 


MOS2 


30 


PPS 


2.0 


SWIRE2 J022255.87-051351.7 


0.846 


1 


0111110501 


20056 


PN 


20 


IXMM 


22 


SWIRE2 J022306.74-050529.1 


0.330 


2 


0111110501 


20056 


PN 


46 


IXMM 


0.95 


SWIRE2 J022348.96-025426.0 


0.355 


Sbt/L 


0037981601 


13319 


PN 


;£ 26 


PPS 


;S 0.8 


SWIRE2 J022356. 49-025431.1 


0.451 


2 


0037981601 


13319 


PN 


i 26 


PPS 


;S 0.8 


SW1RE2 ,7022422.27-031054.7 


1.226 


1 


0037981201 


8420 


PN 


170 


2XMM 


8.4 


SW1RE2 ,7022438.97-042706.3 


0.252 


2 


0112680501 


10040 


PN 


1600 


2XMM 


30 


SW1RE2 ,7022455.77-031153.6 


0.790 


2? 


0037981301 


5702 


PN 


< 10 


PPS 


< 0.8 


SWIRE2 J022542. 02-042441.1 


0.580 


2 


0112681001 


41946 


PN 


;S 442 


PPS 


;S 0.4 


SWIRE2 J022600.01-035954.5 


1.095? 


IR 


0112680201 


13675 


PN 


33 


2XMM 


1.0 


SWIRE2 J022612. 67-040319.4 


0.330 


2 


0112680201 


13675 


PN 


i. 27 


PPS 


;S 0.8 



5.3. Reliability and effectiveness of mid- 
infrared selection of AGN 

The IRAC color technique is clearly effective at 
selecting AGN, at least for objects with fluxes of 
a few mJy at 24/im. Wc have been able to un- 
ambiguously classify most of our objects as AGN 
based on their optical spectra. Contaminants are 
certainly possible, however. At moderate redshift, 
starbursts will be able to move into the AGN 
"wedge" if they have large amounts of warm dust, 
and at high rcdshifts, the 1.6/im bump of stellar 
SEDs can lead to IRAC colors which mimic those 
of AGN. We can state that for this sample of 77 
bright 24/im sources, only seven lack obvious op- 
tical signatures of AGN activity, and of these, one 
is detected in an XMM field, which implies a con- 
tamination rate of < 6/77 (8% ). 

Of course, the selection of bright 24/xm sources 
itself produces a large fraction of AGN, as demon- 
strated by Brand et al. (2006). In particular, when 
combined with optical photometric criteria, a sam- 
ple containing a large fraction of AGN can be 
made by simply rejecting the optically-bright ob- 
jects. To investigate how effective such a cut is, 
and to take a first look at the completeness of our 
AGN sample, we have examined the brightest 50 
sources in the XFLS 24/im list (Table 6). This 



sample contains 29 objects which are classed as 
AGN according to the criteria of Section 3. As the 
starbursts are also concentrated to lower rcdshifts, 
this qualitatively confirms the result of Brand et 
al. (2006), namely that the mid-infrared luminos- 
ity density is dominated by AGN at mJy flux lev- 
els. Objects with faint optical identifications are 
certainly more likely to be AGN - no object with 
i? > 21 is classified as a starburst or LINER, but 
apart from that division there is a wide range in 
magnitudes (i? w 17 — 20) over which AGN and 
starburst/LINERs overlap, particularly if no opti- 
cal morphology criterion is used and thus quasars 
are included. 

This sample of the 50 brightest 24/im sources 
also allows us to take a first look at how complete 
the AGN selection using IRAC colors is. Interest- 
ingly, we find that seven objects outside the IRAC 
color-based AGN selection show AGN-likc emis- 
sion line ratios in a BPT diagram, though these are 
nearly all close to the AGN-starburst transition re- 
gion (Figure 8) . Their positions in the IRAC color- 
color plot are shown in Figure 7. In six of these 
cases IRAC 8.0/im band flux is very high, pushing 
the objects into the starforming galaxy region of 
the color-color plot of Figure 1. The remaining 
object (SSTXFLS J171 150.2-^590041) shows pure 
Seyfert emission line ratios, but has only a weak 
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Table 6: The 50 brightest extragalactic 24;um sources in the XFLS with 4-band IRAC coverage 

Name S24 (mjy) R Spectrum* Redshift Class Basis IRAC-selected AGN? 



SSTXFLS J171115. 2+594907 


9.41 


20.97 


L06 


0.587 


Sbt 


strong Balmer lines 


yes 


SSTXFLS J171150.2+590041 


10.88 


17.97 


P06 


0.060 


2 


BPT 


no 


SDSS J171207.4+584754 


13.34 


17.76 


L06 


0.269 


1 


Fell 


yes 


SSTXFLS J171232. 3+592126 


8.47 


17.56 


P06 


0.21 


Sbt/L 


low ionization 


no 


SSTXFLS J171302. 3+593611 


11.83 


21.50 


L06 


0.668 


2 


[Nev]33 


yes 


SSTXFLS J171313. 9+603146 


10.49 


18.31 


L06 


0.105 


2 


BPT 


yes 


SSTXFLS J171325. 1+590531 


9.43 


18.33 


L06 


0.126 


2 


BPT 


yes 


SSTXFLS J171335.1+584756 


23.66 


17.56 


L06 


0.133 


IR 


Fell 


yes 


SDSS J171352. 4+584201 


23.90 


17.53 


SDSS 


0.521 


1 


BL/Feii 


yes 


SSTXFLS J171414.8+585221 


8.97 


18.70 


P06 


0.167 


Sbt 


BPT 


no 


SSTXFLS J171430. 7+584225 


8.24 


19.72 


L06 


0.561 


IR 


BL 


yes 


SSTXFLS J171437. 4+595647 


11.60 


18.75 


P06 


0.196 


2 


BPT (composite) 


no 


SSTXFLS ,7171446.4+593359 


7.37 


17.67 


P06 


0.129 


Sbl/L 


low ionization 


no 


SSTXLFS ,7171530.7+600216 


11.55 


18.82 


L06 


0.420 


2 


[Ncv]26 


yes 


SDSS J171540. 1+591647 


6.88 


17.43 


SDSS 


0.116 


Sbt/L 


low ionization 


no 


bbTXFLb J171542. 0+591657 


26.33 


17.50 


P06 


0.116 


Sbt 


BPT 


no 


SSTXFLS .7171544.0+600835 


6.86 


17.07 


L06* 


0.157 


2 


BPT 


no 


SDSS J171607. 2+591456 


34.57 


16.27 


SDSS 


0.054 


Sbt 


BPT 


no 


Oi-JOO Jl/lDl'i.4i-t-Oy041ZO 


8 41 


1 7 

1 / .oy 


ouoo 


n 1 K"^ 
u.ioo 


CJKf 


BPX 


no 


Q'RQCI T171fi^n 1_l_fini/10^ 
OiJOiD J 1 ( IDoU. l-|-DU14tZO 


0.00 


1 7 no 




n 1 07 


oDt 


BPX 


no 


00 -L J\.r IjO >J i / iDo4:. U-f-DU 141:^0 


1 n 7fi 


17 04 


rVO 


n 1 ns 


Sbt 


BPX 


no 


SDSS J171635. 9+601436 


7.66 


17.07 


SDSS 


0.106 


Sbt 


BPT 


no 


SDSS J171641.0+591857 


20.88 


16.41 


SDSS 


0.056 


Sbt 


BPT 


no 


iDOlArljiD J 1 / lDOU.O-(-OyO / 01 


D.oy 


1 Si SI 


T OR 


n 1 so 

U. lOZ 


2 


'R'P'T' /'l-<-%TlTr-li-10 1 + £l\ 

1)1 i ^^composixe^ 


yes 


QQT'VPT G T1 71 71 1 1 _LKn071 fl 
00 1 -A.r IjO J±/±/li.i-t-DU.^iiU 


9.43 


17.88 


LOO 


0.110 


QKf /J 
Dt / Li 


low ionization 


no 


GGTVPT Q. TITIT/I/I 1 -L '^S'?S/1 S 

00 1 J\.r LjO >> L I 1 I 41:4:. l-t-00004:0 


1 Q9 

lo.yz 


1 7 7'^ 
±1.10 




066 




BPT 


no 


SSTXFLS ,7171754.6+600913 


8.94 


21.80 


L06 


4.27 


2HZ 


narrow UV lines 


yes 


SSTXFLS ,7171831.7+595317 


8.27 


20.56 


L06 


0.700 


2 


[Nev]26 


yes 


ooiy^rijo ,1 L { LO'i^). i -\-'} J-->'^'} J 


13.22 


17.99 


L06 


0.383 


1 


■RT 


yes 




1 QO 

10. y.^ 


1 Q 9^ 


pnfi 
r^UD 


322 


2/L 


BPX 


no 




Oft Q1 


1 7 1 n 

1 / . lU 


T flR 
ijUO 


n 1 78 




■RT /TTaTT 
JjIj / r GIL 


yes 


SSTXFLS J171913. 5+584509 


8.90 


19.00 


L06 


0.318 


2 


[Nev]12, BPX 


yes 


SSTXFLS J171916.6+593449 


17.40 


17.91 


P06 


0.166 


2/L 


BPX 


no 


SSTXFLS J171933.3+592742 


7.53 


17.71 


POO 


0.139 


2 


BPT (composite) 


no 


SDSS J171944.8+595705t 


6.73 


16.66 


SDSS 


0.069 


Sbt 


BPX 


no 


SDSS J171944. 9+595707+ 


7.44 


16.66 


SDSS 


0.069 


Sbt 


BPX 


no 


SSTXFLS J172028. 9+584749 


9.57 


16.74 


L06* 


0.127 


2 


BPX 


no 


SDSS J172043.2+584026 


9.48 


17.08 


SDSS 


0.125 


Sbt 


BPX 


no 


SSTXFLS J172050.5+591512 


9.41 


21.62 


L06 


1.44 


IR 


BL 


yes 


SSTXFLS J172123. 1+601214 


13.34 


18.78 


L06 


0.325 


2 


[Nev]33, BPX 


yes 


SSTXFLS J172219. 5+594507 


7.76 


18.71 


L06 


0.272 


2 


BPX 


yes 


SSTXFLS J172220. 2+590949 


8.78 


17.58 


P06 


0.179 


Sbt 


BPX 


no 


SSTXFLS ,7172228.0+601525 


7.12 


19.85 


L06 


0.741 


2 


[Nev]15 


yes 


SDSS J172238. 7+585107 


6.73 


18.17 


SDSS 


1.617 


1 


BL 


yes 


SSTXFLS J172328. 4+592947 


8.10 


22.15 


LOO 


1.34? 


IR 


BL 


yes 


SSTXFLS J172400.6+590228 


13.34 


17.99 


P06 


0.178 


Sbt 


BPX 


no 


SDSS J172402.1+600601 


7.76 


17.12 


SDSS 


0.156 


Sbt/L 


low ionization 


no 


SSTXFLS J172508.8+591016 


9.40 


17.78 


P06 


0.164 


Sbt 


BPX 


no 


SSTXFLS J172542.3+595317 


9.45 


20.03 


LOO 


0.437 


Sbt/L 


low ionization 


yes 


SSTXFLS J172619.7+601559 


6.57 


17.13 


LOO 


0.925 


1 


BL 


yes 


■Provenance of spectra: L06 - this 


paper, P06 


- Papovich et al. (2006), 


SDSS 


- Sloan Digital Sky Survey data release 4. 




^ Additional spectra obtained using COSMIC 


on the Palomar 200-inch. 


, 2005, 


June 27-28 







Parts of a single galaxy 
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mid-infrared excess and is in the part of Figure 

1 dominated by quiescent galaxies. This suggests 
that in some cases the signature of AGN activity 
in the mid-infrared can be overwhelmed by strong 
star formation activity, at least ai z ^ 0.3 where 
the PAH features can still contribute to the IRAC 
channel 4 flux. It may also indicate that much star 
formation activity (particularly any close to the 
nucleus) may be more obscured by dust than the 
narrow-line region of an AGN in the same galaxy. 
An AGN selection with the diagonal cut of the 
wedge raised by x +0.5 in lg(S'8.o/'S'5.8) would 
therefore increase the completeness of the AGN 
selection at low redshifts, though at the risk of 
reducing reliability. 

6. The ratio of obscured to unobscured 
AGN 

6.1. Redshift and luminosity distributions 

Figure 9 shows the redshift histograms for our 
various classes of objects, compared to the redshift 
distribution of the sample as a whole. 

Figure 10 shows the corresponding luminosity 
distributions of our objects. We calculate lumi- 
nosities at a rest-wavelength of 5^m. This region 
of the continuum is free from features, and can be 
calculated by interpolation between measurements 
over our full range of redshifts. The most heavily 
reddened AGN will, however, still be significantly 
extincted at this wavelength, which is a possible 
bias in Figure 10. 

As previously discussed, there are strong selec- 
tion effects favoring the selection of low-z, high lu- 
minosity AGN in our sample. Therefore we have 
made no attempt to investigate in detail the ratio 
of obscured to unobscured AGN as a function of 
redshift or luminosity to compare with, for exam- 
ple, the receding torus model of Lawrence (1991). 

Overall, the breakdown into optical spectral 
types is as follows: 33% (25/77) normal type-1, 
14% (11/77) reddened type-1, 44% (34/77) type- 

2 (including 2HZ), and 9% (7/11) starbursts and 
ambiguous starbursts/LINERs. Defining an ob- 
scured quasar as anything with E{B — V) }i 0.3 
in the rest frame and an AGN nature confirmed 
through optical emission lines or X-rays (the type- 
IRs, type-2s, and the single X-ray detected star- 
burst/LINER) we obtain a ratio of obscured to 
unobscured objects of 1.8:1. If we assume that 



Type 1 
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12 3 4 5 
Redshift 



Type 2 



1 2 3 4 5 
Redshift 
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6 - 



1 2 3 4 5 
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Fig. 9. — The redshift distribution for our different 
classes of object. The total distribution for the 
whole sample is shown as the dotted line on all 
the plots, and the solid line the distribution by 
type. 
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Fig. 10. — As Figure 6, but for rest-frame lumi- 
nosity at 5/im. The vertical dashed line is at the 
approximate division between Seyfert and quasar 
luminosities. 
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the other six starburst and starburst/LINER ob- 
jects also contain obscured AGN this increases to 
2.1:1. 

The fraction of type-lRs is higher than that 
found by Brown et al. (2006) in their 24/im se- 
lected sample. Brown et al find 20% of their 
broad-linc objects arc classed as reddened type-Is, 
compared to a type-lR/ (type-l+type-lR) fraction 
of 30% (14/46) in our sample. Our deeper spec- 
troscopy limit {R ~ 23 compared to i? < 21.7) is 
probably responsible for this difference. Two of 
our type-lR quasars have R < 21.7, and a third 
has R = 21.6, excluding these would reduce our 
type-lR fraction to 25%. 

Before we can relate these ratios to an estimate 
of how much supermassive black hole growth is 
obscured, we need to bear in mind the intrinsic 
variations in quasar SEDs. Selection in the mid- 
infrared will inevitably mean our sample is biased 
in favor of objects whose intrinsic SEDs (before 
reddening of the quasar nucleus) are particularly 
bright in the mid-infrared, e.g. because of an un- 
usually large dust covering factor. For type-1 and 
type-lR objects we have good evidence that this 
is not a large effect. Brown et al. (2006) find that 
the space density of their infrared-selected quasars 
is very similar to that of optically-selected sam- 
ples, consistent with the small systematic varia- 
tions seen in the quasar SEDs of Richards et al. 
(2006) and Hatziminaoglou et al. (2005) between 
optically-bright and infrared-bright quasars. Al- 
though our type-1 quasars at z > 1 are redder 
in i?-[8.0] color than the Richards et al. template 
color in Figure 5, most are close to the template, 
and several at z < 1 are bluer. For the more 
heavily obscured objects there are two arguments 
that can be brought to bear. First, even in type- 
1 quasars about 40% of the UV emission is re- 
processed by dust (Richards et al. 2006). Even 
if the obscured quasars have dust covering fac- 
tors near unity, their dust emission would only 
increase by a factor ~ 2 (and they would be signif- 
icantly extincted in the mid-infrared if that were 
the case, leading to their infrared luminosities be- 
ing underestimated). Second, the majority of the 
obscured quasars possess luminous narrow emis- 
sion lines, showing that a significant fraction of 
lines of sight from the nucleus are able to escape 
to ionize the narrow-line gas. Furthermore, in the 
case of radio-loud objects, where the radio lobe lu- 



minosity provides an independent measure of the 
AGN luminosity, Haas et al. (2004b) find that 
radio-loud quasars and their type-2 equivalents, 
radio galaxies, show a statistically-identical ratio 
of radio to 10-1000/zm infrared luminosity, indi- 
cating that the obscured population has a simi- 
lar dust covering factor to the unobscured pop- 
ulation, as expected from orientation-based Uni- 
fied schemes (Antonucci 1993). We are thus con- 
fident that mid-infrared selection docs not signifi- 
cantly bias an AGN sample towards objects with 
unusually high dust covering factors, and so we 
are fairly confident that our AGN sample is a fair 
proxy for one selected on the basis of bolometric 
flux. We can thus safely conclude that ~ 67% 
of black hole growth in the redshift and lumi- 
nosity range probed in this study is obscured by 
E{B -V)> 0.3. 

6.2. Comparison to other Spitzer selected 
Scimples 

Mid-infrared selection of obscured AGN has 

been successfully performed by several other 
groups. Stern et al. (2005) use Hectospec data 
on the Bootes field to study AGN selected using 
a mid-infrared color cut similar to that in Figure 
1, but with different axes. Plotting the objects 
of this paper onto the selection of Stern et al. we 
find that only a few obscured AGN found within 
our mid-infared color criteria are missed by their 
color selection. They flnd an w 3 : 1 ratio of ob- 
scured to unobscured AGN, somewhat higher than 
found in this study, which is 2 : 1. However, as 
Richards et al. (2006) point out, their result relies 
on a power-law interpolation of the quasar SED 
between the mid-infrared and optical, which is 
not quite correct (see, e.g., the variation of R-[8.0] 
color in Figure 5), so their obscured AGN frac- 
tion may be overestimated. Alonso-Herrero et al. 
(2006) use another slightly different mid-infrared 
selection, based on the presence of a power-law 
mid-infrared SED, to find AGN in the Chandra 
Deep Field South. Due to this rather more com- 
plicated selection criterion it is hard to assess 
how many objects in the color-selected AGN sam- 
ples would satisfy the selection criteria of Alonso- 
Herrero et al., but most of them are likely to, as 
our color criterion is very similar to a power-law 
selection. Their study complements the study of 
this paper very well, as their 24/im fluxes range 
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from w 0.05 — 3mJy compared to our range of 
« 4.4 - 20mJy. Their median redshift of ~ 1.4 
(based mostly on photometric redshifts), compares 
with a median of 0.6 for the combined spectro- 
scopic XFLS and SWIRE samples. Their objects 
are thus typically of similar, or only slightly lower 
luminosity to the objects in this study. The fact 
that they find a similar ratio of obscured to un- 
obscured AGN of » 2 : 1 (based on SED fitting) 
is therefore consistent with no strong dependence 
on this ratio on redshift in the range ~ 0.5 — 2. 

Martmez-Sansigre et al. (2005,2006) use a dif- 
ferent approach to selecting their obscured AGN. 
By searching for 24/xm sources with radio emission 
in excess of the radio- infrared correlation they are 
able to isolate obscured AGN irrespective of their 
mid-infrared SED shape. Although the proper- 
ties of these radio-intermediate objects may be dif- 
ferent from the predominately radio-quiet popula- 
tions found from the mid-infrared color selection, 
the ratio of obscured to unobscured objects is sim- 
ilar, in the range 1-3:1. 

In the X-ray, Alonso-Herrero et al. (2006), Don- 
ley et al. (2005) and PoUetta et al. (2006) find a 
significant population (40-60%) of Spitzer selected 
AGN which arc highly-obscured in the X-ray. In 
particular, many 5piteer-selected AGN in the sam- 
ples of Alonso-Herrero et al. and Donley et al. are 
missing from even the ^Ms exposures of the Chan- 
dra Deep Fields. 21 of the objects in our samples 
overlap XMM-Newton fields, we detect 14 of them 
in ^ 10 — 50ks exposures. Given the difference 
in flux limits at 24/im, the relative X-ray depths 
are comparable, although the higher redshifts of 
the Alonso-Herrero ct al. and Donley et al. sam- 
ples should lead to these objects being more de- 
tectable at X-ray wavelengths due to the negative 
fc-correction of an absorbed X-ray spectrum. 

One final point of interest is that the three high- 
est redshift objects in our sample are all obscured 
quasars. Although the statistics are small, the fact 
that these objects were found in spite of strong se- 
lection eff'ects in favor of lower redshift objects is 
an indication that at very high redshifts {z ^ 3) 
obscured accretion may be even more common. 

7. Conclusions 

The results presented in this paper, taken to- 
gether with studies of higher redshift objects se- 



lected from mid-infrared deep fields and objects 
found using joint mid-infrared - radio selection 
show that a large population of luminous obscured 
AGN exist at z > 0.3, including a significant 
population with bolometric luminosities above the 
Seyfert-quasar divide. We find that the obscured 
population dominates in number density over the 
unobscured one over a wide range in AGN lumi- 
nosity, from Seyferts to luminous quasars. Com- 
parison with Alonso-Herrero et al. (2006) suggests 
little variation in the obscured to unobscured ra- 
tio with redshift for AGN with strong mid-infrared 
emission. 

Our results contrast with X-ray studies, which 
typically find that obscured AGN are confined 
to low redshifts (z ^ 1) and luminosities below 
the Seyfert-quasar divide Lx ~ 10'*^ergs~^ (e.g. 
Trump et al. 2006) , suggesting that high luminos- 
ity and high redshift obscured quasars are miss- 
ing from current X-ray based surveys, presumably 
due to limitations of areal coverage and depth. 
Although objects with columns Nh ^ 10^'^cm^ 
can be found in deep X-ray surveys, locally ~ 
50% of AGN are Compton thick, with columns 
> Nh ~ lO^^cm"^, and the equivalents of these 
objects at high redshifts would be very hard to 
find in the X-ray, even allowing for the positive 
k-correction. This has important implications for 
models of the X-ray background (e.g. Worseley et 
al. 2006). However, it is clear that all the tech- 
niques used to find obscured AGN are incomplete. 
Mid-infrared selection presumably fails for weak 
AGN, for highly obscured AGN, and at high red- 
shifts. Our highest redshift objects are selected 
not on the basis of their hot dust emission, but 
on the basis of dust-reddened accretion disk emis- 
sion, so wc will inevitably miss the most obscured 
objects at z > 3. 

For most objects, classification of galaxy type 
(AGN or starburst) based on optical emission 
line ratios and classification based on mid-infrared 
IRAC color appear to agree fairly well. The mid- 
infrared criterion is seen to fail for z ^ 0.3 AGN 
with strong PAH emission in the IRAC 8.0/im 
band from an accompanying starburst. There 
are also some objects with AGN-like IRAC colors 
which lack the high-ionization emission line spec- 
tra of AGN. For most of these, their nature (AGN 
or starburst) remains to be resolved. 

The discovery of a large population of ob- 
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scured AGN has important implications for ac- 
cretion efficiency and black hole demographics. 
The mass density of black holes in the Universe 
today may be accurately estimated using the 
bulge-luminosity black hole mass relation, and 
the amount of accretion luminosity arising from 
unobscurcd accretion may be estimated from the 
optical quasar luminosity function (Soltan 1982). 
Recent estimates using this technique place the ra- 
diative efficcncy of accretion close to that expected 
from a non-spinning (Schwartzschild) black hole 
(e.g. Yu & Tremaine 2002). However, if most ac- 
cretion is actually hidden, to maintain the same 
mass density of black holes today in the face of 
more accretion activity requires a higher radiative 
efficiency, which can only be obtained from rapidly 
spinning Kerr black holes. 
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A. Comparison of classifications with those of Papovich et al. 

In an independent study, Papovich et al. (2006) obtained spectra of bright 24/im XFLS sources with 
the Hectospec Fibre spectrograph on the MMT. Several of these objects are common to our XFLS AGN 
sample. Papovich et al. attempted to classify their objects using a semi-automated algorithm. For the 
most part, our classifications agree, and where we differ it is usually because the Papovich et al. "galaxy" 
classification generally includes narrow- line AGN (Table 6). There are, however, a few objects classed as 
broad line quasars or AGN for which there seems to be no evidence for broad lines in either our spectra or 
the Hectospec spectra. In one case, however, (SSTXFLS J171430. 7+584225) the Hectospec spectrum docs 
show weak broad H/3 and Mgii2798 emission, which appear very weak or absent in our Palomar/COSMIC 
spectrum. The Hectospec spectrum is of higher signal-to-noise than our spectrum, and the H/3 line falls 
close to the 7600A atmospheric absorption feature (for which we do not attempt to correct), so it is most 
likely that we simply missed the weak broad H/3 line in this case. We have therefore reclassified this object 
as type-lR. One further discrepancy is in our classification of SSTSFLS J171454. 4-1-584948, ths object is 
classed at QSO/starburst by Papovich et al., but plots among the starbursts the BPT plot of Figure 4b, and 
has thus been classified as a starburst. 



B. Spectroscopy of the 8/im flux limited sample 

Most of the objects in the 8/xm flux-limited sample of candidate obscured AGN of Lacy et al. (2004, 
2005a) are also members of the XFLS sample of Table 3, however, a few have 24/im fluxes below 4.6mJy and 
are hence missing. For completeness we list their redshifts, where measured, and spectral types in Table 7 
(see also Lacy et al. 2005). 
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Table 7: Non-SDSS objects in the XFLS sample with spectra in Papovich et al. (2006) 



Name 


S'24.o/mJy^ 


Redshift 


Class 


Basis 


P06 class* 


SSTXFLS J171115.2+594906 


9.4 


0.587 


Sbt 




galaxy 


SSTXFLS J171147.4+585839 


4.8 


0.800 


2 


[Nev]13 


galaxy 


SSTXLFS J171324.2+585549 


4.9 


0.609 


2 


[Nev]12 


galaxy 


SSTXFLS J171325.1+590531 


9.5 


0.126 


2 


BPT 


QSO/AGN 


SSTXFLS J171331. 5+585804 


5.9 


0.435 


Sbt/L 




Galaxy 


SSTXFLS J171335. 1+584756 


23.7 


0.133 


IR 


Fell 


QSO/BL 


SSTXFLS J171430. 7+584225 


8.2 


0.561 


2 


[Nev]20 


QSO/BL 


SSTSFLS J171454.4+584948 


4.9 


0.253 


Sbt 


BPT 


Galaxy / starforming 


SSTXLFS J171513.8+594638 


5.0 


0.248 


IR 


BL (radio-loud) 


QSO/starburst 


SSTXFLS J171530.7+600216 


11.6 


0.420 


2 


[Nev]26 


QSO/BL 


SSTXFLS J171650. 6+595752 


6.7 


0.182 


Sbt 


BPT 


galaxy / starburst 


SSTXFLS J171708. 6+591341 


5.4 


0.646 


2 


[Nov] 37 


QSO/BL 


SSTXFLS ,7171839.6+593359 


10.9 


0.383 


1 


BL (radio loud) 


QSO/BL 


SSTXFLS ,7171913.5+584508 


8.9 


0.318 


2 


[Nev]12, BPT 


QSO 


SSTXFLS J172219. 5+594506 


7.8 


0.271 


2 


BPT 


QSO/AGN.BROAD 


SSTXFLS J172228. 1+601526 


7.1 


0.741 


2 


[Nev] 


Galaxy / manual 


SSTXFLS J172245.0+590328 


4.7 


0.797 


2 


[Nev]43 


Galaxy 



* P06 class is the classification of the object in the MMT/Hectospec survey of Papovich et al. (2006), in the form class/subclass. 



Table 8: Objects in the 8/xm selected sample of Lacy et al. 2004 not in the 24/xm flux-limited XFLS sample. 



Name 


Sa.o 
(mJy) 


z 


Class 


Basis 


SSTXFLS J171106. 8+590436 


1.38 


0.462 


2? 


[OinJ5007/[OiiJ3727=3.8 


SSTXFLS J171133.4+594906 


1.91 


7 


1? 




SSTXFLS J171421. 3+602239 


1.25 


? 


2/Sbt/L? 




SSTXFLS J171804. 6+602705 


1.05 


0.43? 


2/Sbt/L? 




SSTXFLS J171930. 9+594751 


1.44 


0.358 


2 


[Nev] 30 


SSTXFLS J172253. 9+582955 


1.00 


1.753* 


1 


BL 


SSTXFLS J172458.3+591545 


1.10 


0.494 


2? 


[Oin] 5007/ [Oii]3727= 1.8 


SSTXFLS J172601.8+601100 


1.43 


1.12? 


IR 


single BL 


spectrum from Papovich et al. 


(2006). 
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